Chimera Ising Walls in Forced Nonlocally Coupled Oscillators 
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Nonlocally coupled oscillator systems can exhibit an exotic spatiotemporal structure called 
chimera, where the system splits into two groups of oscillators with sharp boundaries, one of which 
is phase-locked and the other is phase-randomized. Two examples of the chimera states are known: 
the first one appears in a ring of phase oscillators, and the second one is associated with the two- 
dimensional rotating spiral waves. In this article, we report yet another example of the chimera state 
that is associated with the so-called Ising walls in one-dimensional spatially extended systems, which 
is exhibited by a nonlocally coupled complex Ginzburg-Landau equation with external forcing. 

PACS numbers: 05.45.Xt, 82.40.Ck 



I. INTRODUCTION 

Nonlocally coupled oscillator systems can exhibit a 
remarkable class of patterns called chimera, in which 
identical oscillators separate sharply into two domains, 
one coherent and phase locked, the other incoherent and 
drifting [l|, 0, H, 13] • These two groups of the oscilla- 
tors together maintain stable organized patterns in the 
system. The existence of such patterns was first no- 
ticed and explained in a ring of phase oscillators [l[ , and 
studied in further detail Furthermore, it was found 
that the chimera states also appear in rotating spiral 
waves in two-dimensional spatially extended systems (H • 
Recently, an interesting pattern similar to the chimera 
states studied in Refs. [l|, y| was also found in nonlocally 
coupled Hodgkin-Huxley equations with excitatory and 
inhibitory synaptic coupling, where such patterns spon- 
taneously appears as a result of the instability of the 
uniform state, in contrast to the original chimera that 
appears when the uniform state is stable @. 

In this paper, we present another simple example of 
the chimera state associated with the Ising walls in one- 
dimensional spatially extended systems, which is exhib- 
ited by a nonlocally coupled complex Ginzburg-Landau 
equation (CGLE) with a parametric forcing. Recently, 
a nonlocally coupled CGLE without a forcing has been 
studied intensively [H, ELI 03, GJ, ELI El 51, and a 
locally coupled CGLE with a forcing has also been in- 
vestigated widely [H HH [H [H, 0, [H, [26j]. To our 
knowledge, Battogtokh considered the nonlocally coupled 
CGLE with the forcing for the first time, and demon- 
strated various interesting phenomena, e.g., nonequilib- 
rium Ising-Bloch transitions [6|. Here, we will focus only 
on the chimera Ising wall, which is the simplest one 
among them, but its theoretical analysis has not yet been 
carried out. 

The organization of the present paper is the following. 
In Sec. [Til we introduce our model and illustrate its nor- 
mal and chimera Ising walls by numerical simulations. In 



Sec. IIII1 we reduce our model to the phase model, and 
numerically demonstrate that the reduced phase model 
also exhibits the chimera Ising walls. In Sec. IIV1 a func- 
tional self-consistency equation is derived by introducing 
a space-dependent order parameter, and its numerical 
solution is compared with the numerical simulation pre- 
sented in Sec. IIIII Concluding remarks will be given in 
the final section. 



II. CHIMERA ISING WALLS 

We consider the following equation that describes a 
system of nonlocally coupled limit-cycle oscillators driven 
by a parametric external forcing, 



d t A = (1 + ic ) A-(l + ic 2 ) \A\ 2 A 
+ K (1 + ici) (B - A) +jA*, 



(1) 



which we call a forced nonlocally coupled complex 
Ginzburg-Landau equation. Here cq, c\, C2, K, and 7 are 
real parameters, a complex amplitude A(x, t) represents 
the state of a local limit-cycle oscillator at location x 
and time t, and A* is the complex conjugation of A. The 
quantity B(x, t) represents the nonlocal coupling defined 
by 



B{x,t) = J dx'G{ 



x - x') A(x' ,t) , 



G ( x ) = 2 exp ^ l 2 ^' 



(2) 



(3) 
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where the nonlocal coupling function G(x) is normal- 
ized in the infinite domain. K represents the cou- 
pling strength, and c\ is the phase shift of the cou- 
pling [H, [lj, (H d, M, EI El- The last term 7 A* 
represents the effect of the parametric external forcing 
with almost double the natural frequency, whose inten- 
sity is given by 7, and en — ci stands for the frequency 
mismatch [U|2l|, [2I, [H \M, HI Hi] . In the absence of the 
coupling and the external forcing, i.e., K = 7 = 0, Eq. (JTJ 
is simply given by dtA = (1 + ic^)A — (1 + ic2)|A| 2 A, 
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which is the simplest limit-cycle oscillator called Stuart- 
Landau oscillator [8, 9], so that Eq. (Q]) describes a sys- 
tem of forced nonlocally coupled oscillators. In addition, 
Eq. (TTJ is a normal form that can be derived from a wide 
class of reaction-diffusion systems near the Hopf bifur- 
cation point under particular assumptions by using the 
center-manifold reduction method 

In our numerical simulations, the continuous media of 
size L was discretized using N — 2048 grid points with 
sufficiently small and fixed separation Ax = 0.005, i.e., 
L = NAx, and the Neumann boundary condition (zero- 
flux) was imposed. Our numerical results are unchanged 
if we further increase the number of grid points N or 
the system size L. The initial condition was such that 
A(x < 0,t = 0) = 1.0 and A(x > 0,t = 0) = -1.0. 
In what follows, we fix the parameter values as cq = 
1.0, ci = —0.5, and C2 = 1.0. Note that we set the 
frequency mismatch at zero, i.e., Co — C2 = 0, for the 
sake of simplicity. Furthermore, we fix the ratio of the 
coupling strength to the force intensity as K/^ = 1.5. 

Figure [1] displays a phase portrait (a), a spatial profile 
of the phase (b), and a spatial profile of the modulus (c) 
for the strong coupling case, K = 1.5 and 7 = 1.0, ob- 
tained by a direct numerical simulation of Eq. (TTJ) . The 
phase portrait is given by a set of grid points in the com- 
plex plane each representing the state of a local oscilla- 
tor at a given time. It is found that a normal Ising wall 
appears and all local oscillators are completely phase- 
locked. 

Figure [2] displays a phase portrait (a), a spatial profile 
of the phase (b), and a spatial profile of the modulus (c) 
for the weak coupling case, K = 0.06 and 7 = 0.04. It is 
found that a chimera Ising wall appears, which consists 
of incoherent drifting oscillators near the center (x = 0) 
and coherent phase-locked oscillators in the peripheral 
regions. All the oscillators take almost the full amplitude, 
|A| = 1, so that the system has no phase singularity. At 
the same time, the spatial continuity of the pattern near 
the center (x = 0) is lost, so that a pair of local oscillators 
infinitely close to each other in this region are not always 
close in the state space. 

We can estimate the critical coupling strength below 
which the oscillator at the center of the normal Ising wall 
starts to drift incoherently, leading to the chimera state. 
From the spatial symmetry of the normal Ising wall, the 
values of A and B vanish at the center, i.e., A(x = 0, t) = 
B(x — 0,t) = 0. If we regard B as an external forcing 
the linearized equation for the complex amplitude 
A(x = 0, i) at the center is given by 



d t A = (1 + ic a ) A - K (1 + ici) A + 7A* 



(4) 



A linear stability analysis of the stationary solution 
A(x = 0, t) = gives the following eigenvalues: 



A± = l-K±J 1 2-(c -Kc 1 y 



(5) 



Therefore, the necessary condition for the appearance of 
the chimera Ising wall, i.e., the condition for the oscilla- 
tor at the center (x = 0) to drift, is expressed as A+ > 0. 
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w> 
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FIG. 1: Phase portrait (a), spatial phase profile (b), and spa- 
tial modulus profile (c) of the normal Ising wall exhibited 
by the forced nonlocally coupled complex Ginzburg-Landau 
equation fTJ in the strong coupling case, K = 1.5 and 7 = 1.0. 
The numerical data are represented by the open circles (only 1 
in every 8 oscillators is plotted). The limit-cycle orbit \A\ — 1 
of the local oscillator is also displayed in (a). 



In fact, when K < 1 is satisfied under our parameter 
conditions, the oscillator at the center (x = 0) starts to 
drift and a chimera Ising wall appears. Hereafter, we fo- 
cus on the chimera Ising walls and consider the situation 
where the coupling strength is much smaller than 1, i.e., 
#<1. 
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FIG. 2: Phase portrait (a), spatial phase profile (b), and spa- 
tial modulus profile (c) of the chimera Ising wall exhibited 
by the forced nonlocally coupled complex Ginzburg-Landau 
equation (JTJ in the the weak coupling case, K — 0.06 and 
7 = 0.04. The numerical data are represented by the open 
circles (only 1 in every 8 oscillators is plotted). The limit- 
cycle orbit \A\ = 1 of the local oscillator is also shown in 
(a). 



III. PHASE REDUCTION 

In order to investigate the nature and the origin of the 
chimera Ising walls in further detail, we reduce Eq. (JTJ) 
to a phase equation for the phase 4>(x, t), which is much 



easier to analyze. When the coupling strength K and 
the forcing intensity 7 are sufficiently small, the phase 
reduction method is applicable, which is actually 

the case for our weak coupling case, K = 0.06 and 7 = 
0.04. The reduced equation takes the form 



d t 4> = -Ka I dx' G(x- x') 
-7Tsin(20 + /3), 



sin 



b' + a)- 



S111Q 



(6) 



where the condition cq — ci — was used. Here 4>' is the 
abbreviation of </>(a/,t). The new parameters are related 
to the original parameters through 



a exp (ia) = (1 — ic\)(\ + 1C2), 
rexp (i/3) = (1 + ic 2 ), 



(7) 
(8) 



where a, r, a, and (3 are real. When the time scale is 
changed so that the force intensity is normalized, t — > 
t/7T, Eq. ([5]) can be rewritten as 



dt<p = — K~ ( J dx' G (x — x') sin( 
-sin(20 + /3), 

where 



sin a 



K~, 



Kg 



(9) 



(10) 



Figure [31(a) displays a spatial profile of the phase (f) 
of the local oscillators obtained by a numerical simula- 
tion of Eq. ([9]) using parameter conditions corresponding 
to Fig. [21 The phase pattern is very similar to that in 
Fig. [2b) . As in the previous case, only the local oscil- 
lators near the center (x = 0) are drifting incoherently, 
and all other oscillators are phase-locked. 

Now let us introduce a space-dependent complex order 
parameter with modulus R(x) and phase $(x) through 

R (x) exp(i<!> (x)) = dx' G (x — x') exp(i0 {x' , t)) . 



(11) 

We assume the order parameter to be time independent, 
which will be confirmed below. In terms of this order 
parameter, Eq. ^ may be expressed in the form of a 
single-oscillator equation 



d t 



i?sin( 



-$ + a) — sin a 



sin(2(/> + /?), (12) 



or, if we further introduce a space-dependent effective 
force function H x ((f) through 



H x (0) = A% 
Eq. 



i?sin( 



$+aj —sin a 

can be expressed as 

9t4> = -H x (0) . 



-sin 



(20+/J), (13) 



(14) 
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Our effective force function is composed of the internal 
force (the first harmonic term) and the external force 
(the second harmonic term), which is similar to the order 
function introduced by Daido [l(| HH, 03 • 

The spatial profiles of R(x) and Q(x) obtained by the 
numerical simulation of Eq. ([9|) are displayed in Fig. [31(b) 
and Fig. (3£c) , respectively. We can see that the order- 
parameter modulus has a vanishing value at the center, 
and both the modulus and the phase of the order param- 
eter are time independent. The distribution of the mean 
frequency Q{x) of the local oscillators, which is defined by 
a long-time average of d t (j), is also displayed in Fig.[3fd). 
In the phase-locked domain, the oscillation frequencies 
are identically zero, while in the drifting domain they 
are distributed. 

The dependence of the effective force function H x ((j>) 
on x and 4> is displayed in Fig. [4] It is found that H x ((j>) 
crosses the zero plane only twice in every 2ir interval of <p 
at each x in the peripheral regions, which will make our 
analysis simple. We should also note that H x ((f>) is always 
below zero near the center {x — 0), namely, H x ((f>) never 
crosses the zero plane in this region. In the next section, 
we develop a self-consistent theory for determining the 
spatial profiles of the quantities displayed in Fig. [3l 



IV. SELF-CONSISTENT THEORY 



(a) % 









A' 



(b) 




We now develop a self-consistent theory that can re- 
produce our simulation results by g eneralizing the ear- 
lier theories [H, 0, i, i, 0, 1 !, [Hi GH, • There are 
two possible cases regarding the solution of Eq. (fl4|) . Let 
-H r min(^) and -ff max (x) denote the minimum and the max- 
imum of H x {(j)) in every 2tt interval of <fi at each x, re- 
spectively. Then, the two cases are expressed by 



H min (x) < < ff max (x) (case I), 



(15) 



and 



Hmi n (x) > or 



0>H max (x) (case II). (16) 



Correspondingly, the oscillators are divided into two 
groups. In the case I, which corresponds to the group 
of the phase-locked oscillators in the peripheral regions, 
Eq. (fl~4"|) admits only one pair of stable and unstable fixed 
points. We denote the stable fixed point formally as 



0o (*) = H- 1 (0) . 



(17) 



where H x x is the inverse function of H x . The average fre- 
quency d){x) of the oscillators in this group is identically 
zero, 



U! {x) = 0. 



(18) 



The case II corresponds to the group of the drifting 
oscillators, for which Eq. (fT4")) admits a drifting solution. 
The average frequency ui{x) is formally expressed as 



u) (x) = 2ir 



2tt 



d(t> 



-H x (0) 



(19) 



(c) n 



O 
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FIG. 3: Spatial profiles of the local oscillator phase <j> (a), the 
order-parameter modulus R (b), the order-parameter phase 
$ (c), and the mean frequency lj of the local oscillators (d), 
obtained by a numerical simulation of the phase equation Q . 
The numerical data are represented by open circles (only 1 in 
every 8 oscillators is plotted). 
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obtained numerically by an iteration procedure, which is 
compared with the results of the numerical simulation 
in Fig. 03 The agreement between the theory and the 
numerical simulation is excellent for all quantities. 

We finally make a brief comment. The order- 
parameter phase does not drift in our chimera Ising 
walls, which makes our analysis easy because we do not 
need to calculate the collective frequency. In the two 
early examples of the chimera states, we had to work 
with a nonlinear eigenvalue problem, or to determine the 
various quantities together with the collective frequency 
(i.e., the eigenvalue) [3, 0. 



FIG. 4: Dependence of the effective force function H x (4>) on 
x and <j>. The contour lines corresponding to H x (<f>) = are 
also displayed in the base plane. Note the gap near the center 

(x = 0), where H x (cj>) does not cross the zero plane. V. CONCLUDING REMARKS 



which depends on x. The contribution to the order pa- 
rameter from the drifting oscillators can be computed by 
the standard method [TJ, @]. That is, we use the invari- 
ant measure, i.e., the probability density p x {4>) associ- 
ated with the drift motion. Noting that the probability 
density for the phase at x to take on value </> must be 
inversely proportional to the drift velocity given by the 
right hand side of Eq. lfT4"|) , we obtain 

=C X [-H x (cj>)]- 1 , (20) 
where C x is the normalization constant given by 



-H x (0) 



(21) 



Putting together the two types of contributions to 
the order parameter, we finally obtain a functional self- 
consistency equation in the form 

S (,M» W )-/tf< I (.-^ W , (22) 
where 

(easel), 

h(x) = { f 2 * .. (23) 

' ' Px ((j))e^d(l) (case II). 



We present a new type of the chimera states associated 
with the Ising walls in the one-dimensional forced non- 
locally coupled complex Ginzburg-Landau equation. We 
have focused on the weak coupling case, where the phase 
reduction method is applicable, and reduced the origi- 
nal model to the phase model. Generalizing the previous 
theories, we derived a functional self-consistency equa- 
tion to be satisfied by the order parameter. Its solution 
successfully reproduced our simulation results carried out 
on this phase model. 

The chimera state studied in Refs. [J Q seems rather 
special because the boundary effects are crucial Q. In 
contrast, the chimera Ising walls stably exist regardless 
of the boundary effects, and they survive even in spa- 
tially infinitely extended systems, like the chimera spi- 
ral waves that persist in two-dimensional spatially ex- 
tended systems [H, Q ■ Our preliminary analysis suggests 
that a chimera hole solution does not seem to exist in 
a one-dimensional nonlocally coupled (cubic) complex 
Ginzburg-Landau equation without a parametric forcing. 
Therefore, our chimera Ising wall may be the simplest ex- 
ample of the chimera states which is still relevant to the 
real-world phenomena. 
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FIG. 5: (Color online) Comparison between the theory and 
the simulation results. The open circles and the solid lines are 
the simulation data and the theoretical curves, respectively. 
Spatial profiles of the local oscillator phase <j> (a), the order- 
parameter modulus R (b), the order-parameter phase $ (c), 
and the mean frequency u) of the local oscillators (d). 



